Spin-charge-density wave in a squircle-like Fermi surface for ultracold atoms 
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We derive and discuss an experimentally realistic model describing ultracold atoms in an optical 
lattice including a commensurate, but staggered, Zeeman field. The resulting band structure is quite 
exotic; fermions in the third band have an unusual rounded picture- frame Fermi surface (essentially 
two concentric squircles), leading to imperfect nesting. We develop a generalized SO(3, 1) x SO(3, 1) 
theory describing the spin and charge degrees of freedom simultaneously, and show that the system 
can develop a coupled spin-charge-density wave order. This ordering is absent in studies of the 
Hubbard model that treat spin and charge density separately. 



Introduction Ultracold atoms in optical lattices have 
recently emerged as a class of condensed matter systems, 
where the properties of the many-body Hamiltonian are 
under exquisite experimental control. Interfering laser 
beams in one, two or three dimensions (D) create stand- 
ing waves: nearly perfect optical lattices for atoms with 
lattice spacing and topology set by the laser geometry 
and wavelength pQ. Optical lattices not only allow for 
the implementation of different lattice models without 
defects, but also open a wide range of possibilities to ma- 
nipulate the parameters of the model describing ultracold 
bosons, fermions, or mixtures thereof. For example, the 
hopping parameters, local chemical potential, and often 
even the interaction strength can be tuned at will. 

Most optical lattice experiments use atoms in a sin- 
gle state [2], however, some experiments study mixtures 
of atoms in two or more atomic "spin" states, each of 
which can experience different lattice potentials [3H5]. 
We derive a lattice model, equally applicable to bosons 
and fermions, with an effective Zeeman magnetic field 
including a term alternating in sign on a site-by-site ba- 
sis |B]. In condensed matter systems, the Zeeman field 
couples strongly to electrons near the Fermi surface [7], 
and in more orchidaceous situations, it breaks local time- 
reversal invariance in topological insulators [8, 9 . 

For particles with two spin states, our lattice model has 
four low-energy bands, and the third is shaped as a squar- 
ish, deformed, Mexican hat for a wide range of system pa- 
rameters. By filling the system with fermions, we obtain 
a peculiar Fermi surface, consisting of the boundaries of 
a squarish ring, essentially two concentric squircles [TO] . 
The particular shape of the Fermi surface suggests that 
nesting effects should be expected. To account for inter- 
actions, we develop an SO (3, 1) x SO (3, 1) description of 
the charge and spin degrees of freedom. Imperfect nest- 
ing along the diagonal connecting corners of the Fermi- 
surface gives rise to a coupled spin-charge-density wave 
(SCDW) instability at a critical interaction strength U c . 
We calculate the imaginary part of the trace of the ran- 
dom phase approximation (RPA) susceptibility to study 
the collective excitations of the system. At the interac- 
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FIG. 1. (a) Schematic layout: two nearly degenerate counter- 
propagating lasers differ in frequency by a very small Sou, and 
are linearly polarized in the x-y plane. An rf magnetic field 
B r { cos(6ujt + cj)) is polarized along z. A bias field Box brings 
the two Zeeman-resolved vtlf levels nearly into resonance, (b) 
Illustrative level diagram: the two wif states are coupled by 
optical and rf magnetic fields. The lattice potential formed by 
the retro-reflected lasers is omitted, (c) Band structure com- 
puted for Qr = 2£, Q r f = 4£, and <j> — 7r/4: the concentric- 
squircles of the Fermi-surface are obtained by filling the sys- 
tem to the third band, and are depicted by the white contours 
at the Fermi energy (also shown is the projection onto the &+- 
k- plane). 

tion strength U C1 a soft mode arises at the optimal nest- 
ing wave- vector Q. The SCDW instability is in general 
incommensurate with the lattice, and is tunable by ex- 
ternal parameters. In contrast with the usual behavior in 
ID, our results show that in 2D a combined treatment of 
spin and charge degrees of freedom is essential to capture 
the possible instabilities of the system. 

The physical system under study [Figs. u\ a) -(b)] con- 
sists of a sample of ultra cold atoms illuminated by two 
pairs of counter-propagating lasers with angular frequen- 
cies uo and uo + Suo (where uo ^> Suo); a third pair of lasers, 
not shown, propagate along ±z and create a ID lattice, 



confining motion to the x — y plane. Our model includes 
a static magnetic field Bq along x and an rf magnetic 
field B v { with angular frequency 5uo along z; rf coupling 
between atoms in spin dependent lattices has been stud- 
ied both experimentally [5] and theoretically [11], where 
the resulting non-trivial real-space lattices suggested po- 
tential application to many body systems and quantum 
computation. In our case, the spin dependence results 
from the interplay of the laser and rf-magnetic fields. 

As was observed in Refs. [6j [12j [13], conventional 
spin independent (scalar) optical lattice potentials ac- 
quire additional spin-dependent terms near atomic res- 
onance: the rank-1 and rank-2 tensor light shifts [15] . 
In the case of alkali atoms, adiabatic elimination of 
the angular momentum J = 1/2 (Dl) and J = 3/2 
(D2) excited states yields an effective Hamiltonian Hq = 
U s (e* -e) + iU v J-(e* xe)/ftforthe J = 1/2 ground state 
atoms, e is the polarization vector of the optical electric 
field and the magnitude of the scalar and vector light 
shifts are related by U v = — 2L 7 " s Aps/3(cj — ojq). Here, 
the fine-structure splitting is Aps — ^3/2 — ^1/2 5 ^1/2 
and fkj 3 /2 are the Dl and D2 transition energies; and 
(Jo = (2cji/2 +^3/2)/3 is their suitable average. U v and 
U s can be independently specified with informed choices 
of laser frequency 00 and intensity. We focus on a prac- 
tical case, where the lasers are detuned far below atomic 
resonance ujo — uo ^> Aps, minimizing spontaneous emis- 
sion and implying \U S \ ^> \U V \ and U s < 0. We express 
momentum and energy in dimensions of Hk r = huo/c and 
E r = h 2 k r 2 /2m, the single-photon recoil momentum and 
energy, respectively, with m the atomic mass. 

The atomic Hamiltonian for the laser and magnetic 
fields in Figs. [lta)-(b) is Ho = U s (cos 2 k r x + cos 2 k r y) + 
9^bJ ' B e ff with B e ff = Box + B r f cos(5ut + <p)z + 
z(U v /2gfjLB) cos(dujt) cos(k r x) cos(k r y) : where the vector 
light shift acts as an effective magnetic field and Bq ^> 
B Y £,U v /2gnB- Here, /ie is the Bohr magneton and g is 
the Lande ^-factor. We select x as the quantizing axes, 
transform into the frame rotating at 5ui, and make the 
rotating wave approximation to find 
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B e ff • z reaches its extrema on the sites of the optical lat- 
tice, giving a bias plus staggered Zeeman field. This pro- 
posal requires the simple retro-reflection of the existing 
"Raman" lasers discussed in Ref. [14 , which were used 
to create an artificial magnetic field (there, B T { was used 
only for state preparation). When \U S \ ^> \U V \ , |<7/iei?rf | 
the conventional tight-binding model [15 , valid when 
U 3 > 5E r , is slightly modified by the effective magnetic 



field evaluated on the lattice sites, yielding 
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Cj jS is an annihilation operator (bosonic or fermionic) 
on site j with spin s; the hopping matrix element t 
can be computed from the band structure of a sinu- 
soidal lattice (for a U s = 5E r scalar lattice t ~ 0.07 'E r )\ 
^rf = g^B^xi\ and £7r = U v /2. Since we focus on very 
small Qr ~ t, the detuning from atomic resonance can 
be quite large. For 40 K, with A F s/2tt = 1.7 THz, the 
detuning is (ujo — oj)/2tt « 50 Aps = 86 THz, yield- 
ing a laser wavelength « 980 nm, far detuned from the 
770.1 nm (Dl) and 766.7 nm (D2) transitions. 

To determine the single-particle spectrum, we define 
spinor field operators Cj = (cj,t? c j,|) T and three com- 
ponent vectors Sj = (S?,Sy,S?) T = cj<7Cj/2, where 
a is the vector of Pauli matrices. Owing to the stag- 
gered Zeeman field, we introduce sublattices A+ and A_ , 
where A± = {(j x ,j y )\(—l) jx + j y = ±1} and we define 
a^± = Cj for j G A±. In addition, we introduce vec- 
tors B± = (O r fcos(0) ± £1r, — Q r f sin(0),O) T describing 
Zeeman fields on the A± sublattices. In this notation, 
the bare Hamiltonian is 

ff = -tX)( a t+ a J.-+ h - c -)+ E s f B + E S J' B - 

In terms of momentum field operators ^ = 
(a£ _j_ ^, a^ _j_ 1 , a^ _ ^, a^ _ , ), the Hamiltonian becomes 
#0 = E k ^]L H o k^k where 



H 



0k 



a • B+/2 -t 7k I 
-t7 k I <r-B_/2 



Here, 7k = 4cos(/c + /2) cos(fc_/2), with k± = ir(k x ± 
k y )/k r , k = (fc_|_,fc_) and I is the 2x2 identity ma- 
trix. The summation goes over the entire Brillouin zone 
— 7r < fc + ,fc_ < 7r, i.e., 7k > 0. The four eigenval- 
ues of H ok are 4 = (t 7k ) 2 + (^ r f/2) 2 + (^r/2) 2 ± 
^rf\/(^7fc) 2 + (^R cos (^)/2) 2 ; together these eigenvalues 
constitute four bands [Fig. (Tic)]: the lowest band has a 
minimum at k = 0, whereas the third band can be shaped 
as a squarish, deformed Mexican hat. The second band 
may either exhibit the same trivial behavior as the first 
band, with a global minimum at k = or have lines of 
degenerate minima along a square contour at the edge of 
the square Brillouin Zone. The fourth band always has 
lines of minima at the zone boundary. 

Spin-charge-density-wave While the single-particle 
spectrum is valid for fermions and bosons, we now fo- 
cus on spin- 1/2 fermions with a Hamiltonian 



H = H + H int ; H ivt = UJ2 c ],t c L c J. 
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FIG. 2. (a) Largest eigenvalue of the static susceptibility ?7Xk(0) as a function of fc+ and k-. The peaks indicate the Q vectors 
for SCDW instabilities, and the largest peak marks the location of the most prominent nesting vector. Imaginary part of 
the trace of the susceptibility Tr[Imxk(^)] in logarithmic scale in the k y -uj plane, with k x = 0, i.e., /c + =— k-. (b) Without 
interactions a linearly dispersing sound mode is observed for small k y . (c) For U/t = 2.60(3), spectral weight builds up for a 



second linear-dispersing mode, which starts from uo — at nk y /k r — 7r/4. (d) For U/t 
intensity at uj = with 7vk y /k r = 7r/4 signals the onset of SCDW instability. 



2.89(3), the sharp increase of the 



The interaction strength U is proportional to the s-wave 
scattering length, and the Fermi energy is chosen to be in 
the third band. The resulting squarish Fermi surface is 
depicted by the white contours in Fig. [lie) and is shaped 
like two concentric squircles. Nesting and local fermion- 
fermion interactions lead to spin- and charge- ordered 
phases in this system. We anticipate a second order phase 
transition when the coefficient of the second order term 
in the Landau free energy vanishes. Formally, we use the 
Hubbard- Str at onovich transformation to treat the inter- 
actions within a saddle point approximation (analogous 
to the time-dependent Hartree-Fock approximation). 

Coexisting spin- and charge-density waves (SDW and 
CDW) have been studied in quasi-2D organics (a few 
coupled chains) using a Monte Carlo aproach [16]. Con- 
ventional approaches to study SDW instabilities in the 
2D Hubbard model neglect the contribution of charge 
density fluctuations [13 [18], which are important here. 
In the following, we develop a generalized solution of 2D 
tight-binding models with local interactions and obtain 
a theory of SCDW instabilities. 

In the coherent states formalism, the grand-canonical 
partition function is Z = J d[c^]d[c]e~ s ^ c ^l h , where 



Hq + H n 



is the Eu- 



We express the inter- 

4 1 c J,i c J,t = 
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S = /c , dr[£ j cj(/a-/i)c j + 

clidean action and j3 = l/fc^T. 

action term in a SO (3, 1) invariant form d 

(l/8)n?-(l/2)Sj T -Sj, with n$ = c]q. The SO(3) imparl- 
ance is required by rotational symmetry; the fact that 
the spin and charge terms have different signs reflects 
the Pauli principle, which requires a vanishing self-energy 
for a polarized state. The Hubbard- Str at onovich trans- 
formation renders the action quadratic in the fermion 
operators by introducing auxiliary bosonic fields, p± 
and M±, which couple to charge and spin density, re- 
spectively. For repulsive interactions the charge den- 
sity term leads to a divergent integral. We resolve 
this problem by integrating along a contour parallel 



to the imaginary axis for p±. Next, we introduce a 
source field J that couples to the charge and spin den- 
sities at each sublattice, and an eight-component vector 
M k , n = (p k5+5n ,M k5+5n ,p k _ 5n ,M k _ 5n ) T expressed in 
terms of momentum k and Matsubara frequency uj n = 
7r(2n + l)/hfi. After integrating out the fermionic fields, 
we obtain a path-integral over the auxiliary bosonic 
field M k?n , which we evaluate in the saddle-point ap- 
proximation. Notice that the saddle-point (M k?n )j de- 
pends on the source field J. We find that d(M)j = 
HU Xj PA dJ, where the generalized RPA susceptibility 
(Xj PA ) _1 = Xj 1 — VU is an 8 x 8 matrix. The ma- 
trix r] = Diag(— 1, 1, 1, 1, —1, 1, 1, 1) is a metric signature 
corresponding to the SO(3, 1) x SO(3, 1) group and xj 
is the bare susceptibility for the renormalized Hamilto- 
nian theory. By neglecting second order fluctuations of 
the Hubbard- Str at onovich fields, we obtain dln(Z[J]) = 
Z7 -1 (M)j • dJ. Neglecting now third- and higher-order 
terms in J we find Z[J] = Z[0] exp(— <S e ff[J]/ft), where 
<S eff = (-h/U)(M)l ■ J - (ft 2 /2)J t • (I - UxoV^Xo ■ J- 
The free energy may be determined by performing a Leg- 
endre transformation f3F[(M)j] = /7- 1 (M)^J-ln(Z[J]) 
(see Ref. [18]) which, up to quadratic order in the devi- 
ation A(M)j = (M)j — (M)o and without an additive 
constant reads 



/^[(M). 



2hU 2 



(A(M)jy ■ ( X « 



RPA\ 



A(M) J . (2) 



The susceptibility x? PA I s evaluated in the absence of 
the source field, J = 0. For homogeneous phases, the 
susceptibility and the Hamiltonian become diagonal in 
momentum and frequency space. Thus, the Hamiltonian 
in the saddle-point approximation becomes 

H k = H 0k -^- J2 P r ^'Tr[P'-'u q l( a )ut 1 ]n F (4 a )). 



q,r,r 7 ,Q; 



(3) 
[Diag(I,0), Diag(<r*,0), Diag^.O), Diag(<j z ,0), 



Diag(0,I), Diag(0, a x ), Diag(0, c^), Diag(0, a z )] T are 
constant 4x4 matrices; N is the number of sites in a 
sublattice; tif(z) = (eP z + 1) _1 is the Fermi distribution 
function; energy is measured with respect to the chemi- 



cal potential e^ — 
which diagonalizes 



.(«) 



/i, and U q is a unitary matrix 



H k = ^UkI (a) uie : 



k ' 



(4) 



with I M = Diag(^ a? i,^ a5 2,^ a ,3,^a,4)- Solving Eqs. (|3| 
and Q self consistently, we determine the RPA suscep- 
tibility (x? PA ) _1 = Xo" 1 — VU, where the k— dependent 
susceptibility at zero source xo is 



Xk ,r (ihfi n ) 



E 



p,a,/3 ' 



^F( £p+k)- n ^(4 /3) ) T r,r / ;q,/3 
P+k,p ' 



£ p+k 



.?(«_ 



ihQ n 



(5) 



*£#:*/ = ^Tr[P'-U p+k I ( ^U^ +k P'-'U p l(«Ut ) ]. 

Here Q n = 2irn/hp is the bosonic Matsubara frequency. 
Static susceptibility The expression in Eq. ([5| can be 
evaluated numerically with arbitrary precision. Though 
our approach is applicable for any temperature regime, 
we restrict ourselves to temperatures close to zero. First, 
we consider the static susceptibility Xk(0) and look for 
possible instabilities. The instability condition for re- 
pulsive interactions requires U > Z7 C , where U c is de- 
termined by det(x^" 1 — U c rj) = 0. Since we avoid the 
van Hove singularity, the susceptibility is finite and the 
critical value U c is nonzero. It is related to the largest 
eigenvalue Aq of the matrix r?Xk(0) by U~ x = maxk Ak- 
Thus, the instability condition becomes Aq/7 > 1, anal- 
ogous to the Stoner criterium. Fig. p[a) shows Aq for 
the Fermi surface in Fig. [I] (0 = 7r/4, Q r = 2t, Q r f = 4t, 
/i = t, and k^T = 10~ 3 t). We calculated Xk(0) numer- 
ically on each point of a mesh with 240 x 240 points. 
The peak with Aq = 0.345(3), corresponding to the crit- 
ical value of interactions U c /t = 2.89(3), is located at 
Q = (7r/4,7r/4), where we expect an imperfect nest- 
ing between inner and outer lines of the Fermi surface 



For these system parame- 



[Fig.|I|:c)] : with4 3 > Q «4 3 >. 

ters, tne eigenvector Vq corresponding to this eigenvalue 
has an anti-ferrimagnetic character and is a mixture of 
both SDW and CDW, hence a SCDW. The details of the 
mixture are not universal. 

The period of the SCDW 2tt/|Q| is in general incom- 
mensurate with the lattice period and is freely tunable by 
changing the vectors B±, and the chemical potential /a. 
Nesting can also occur for other momenta, which results 
into smaller peaks forming the pattern shown in Fig.[2|a). 

Had we neglected the coupling with charge and con- 
sidered only the spin susceptibility, we would find at 
the same value of Q a much lower value for the criti- 
cal interaction strength: U c /t = 2.13(3) compared with 
U c /t = 2.89(3) in the full calculation. In addition, when 
considering only charge excitations, no CDW instability 



occurs for repulsive interactions U > 0. Thus, the cou- 
pling of charge and spin excitations, as developed here, is 
essential to the realization of a phenomenon which other- 
wise would only occur for attractive interactions U < 0. 

Collective excitations Equation ([5| allows us to study 
the collective excitation spectra by analytically continu- 
ing iVt n — )> uo + in and looking at the imaginary part of 
the trace of the RPA susceptibility Tr[Imx k lPA (co')]. For 
k = 10 -2 we find a linear dispersion spectrum in Fig.^Tb) 
in the long wavelength region (the Landau zero sound), 
which could have been anticipated, since we are consider- 
ing a compressible zero-temperature Fermi liquid. At the 
interaction value U c /t = 2.89(3) a soft linearly dispersing 
mode starting from k = Q appears, signaling the onset 
of instability [Fig. |2[d)]. In 40 K, the collective excitation 
spectrum can be experimentally studied with an atomic 
analog of angle resolved photoemission spectroscopy [19] 
or by measuring the dynamic structure factor with en- 
ergy and momentum sensitive Bragg spectroscopy [20] . 

Conclusions We showed how to construct a system 
with a unique Fermi surface consisting of concentric 
squircles. The system has peculiar collective excitations, 
which we analyze in the RPA including both charge- and 
spin- density excitations. Our studies predict an anti- 
ferrimagnetic-like instability combining both CDW and 
SDW with a tunable incommensurate wave- vector, deter- 
mined by the nesting properties of the Fermi surface, for 
sufficiently strong interactions. Moreover, we find that 
the usual approach - neglecting the coupling with den- 
sity fluctuations - significantly underestimates the criti- 
cal value of the interaction strength. 
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